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Abstract 



A generalization of the semiclassical approximation is suggested allowing 
for an essential extension of its region of applicability. In particular, it be- 
comes possible to describe Bose-Einstein condensation of a trapped gas in 
low-dimensional traps and in traps of low confining dimensions, for which 
the standard semiclassical approximation is not applicable. The results of the 
modified approach are shown to coincide with purely quantum-mechanical cal- 
culations for harmonic traps, including the one-dimensional harmonic trap. 
The advantage of the semiclassical approximation is in its simplicity and gen- 
erality. Power-law potentials of arbitrary powers are considered. Effective 
thermodynamic limit is defined for any confining dimension. The behaviour 
of the specific heat, isothermal compressibility, and density fluctuations is 
analyzed, with an emphasis on low confining dimensions, where the usual 
semiclassical method fails. The peculiarities of the thermodynamic character- 
istics in the effective thermodynamic limit are discussed. 
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I. INTRODUCTION 



Physics of Bose gases, exhibiting Bose-Einstein condensation, is currently a subject of 
intensive investigations, both experimental and theoretical (see, e.g., the review works [1- 
5]). A very efficient method for describing Bose-Einstein condensation of trapped atoms 
is the semiclassical approximation, which has been employed for considering the properties 
of Bose-condensed ideal gases trapped in power-law potentials [2,6-11]. The advantage of 
the semiclassical method is its simplicity, as compared to quantum-mechanical calculations, 
and its generality allowing for the treatment of power-law potentials of arbitrary powers, 
except those cases, when the effective confining dimension becomes low. The latter happens 
for the gas of low spatial dimensionality and when the power of the confining potential is 
large. For instance, the semiclassical approximation predicts the absence of Bose-Einstein 
condensation at finite temperature in a one-dimensional harmonic trap [2,7,11], while this 
exists in a quantum-mechanical picture [12]. 

In the present paper, a modification of the semiclassical approximation is advanced, 
which makes it possible to essentially extend the region of its applicability. Thus, the 
low-dimensional gases can be successfully described, with confining potentials of arbitrary 
powers. The definition of an effective thermodynamic limit for trapped atoms, known for 
the case of a three-dimensional harmonic trap [3], can be generalized to arbitrary space 
dimensionality and any confining powers. The specific-heat discontinuity at the condensation 
point can now be studied for all those cases, for which the standard approach does not work. 
The behaviour of isothermal compressibility and the peculiarity of density fluctuations of a 
finite number of trapped atoms, not studied earlier, are investigated. 

Throughout the paper, the system of units is employed, where the Planck constant H = 1 
and the Boltzmann constant k B = 1 are set to unity. 



II. SEMICLASSICAL DENSITY OF STATES 

Let us start with very briefly recalling the basic notions and introducing the necessary 
notation that will be used in the following sections. We shall consider the ideal Bose gas 
confined by means of a trapping potential U(r). The Cartesian vector r = {r a }, with 
a = l,2,...,d, is defined in a rf-dimensional space, so that d = 1,2,.... The trapping 
potential is assumed to be slowly varying in space, such that its characteristic length l be 
much larger than the thermal wavelength X T , 




(1) 

where m is atomic mass and T, temperature. Condition (1) can be rewritten as an inequality 
for the characteristic frequency of the trapping potential uq, 



« 1 . = — 2 ■ (2) 



Under condition (1) or (2), the semiclassical approach is applicable [2,6-10], described 
by the density of states 
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where d > 1 and the integration is over the available volume 

Y £ = {r| U(r) < s} . 
In particular, for a one-dimensional system, one has 

27r yJe-U(x) 

with the turning points ±x e given by the equality U{x e ) = e. Note that Eq. (4) follows 
exactly from Eq. (3). 

The general form of a power-law confining potential can be represented as 



d/2 



a=l 



UJ a 






2 


l a 





L = -== , (5) 



where (v a ,l a , and n a are positive parameters. The characteristic trapping frequency and 
potential length are given by the geometric averages 

/ d \V<* / d \V<i 

= n ^ J ' ^o=^n^J • ( e ) 

An important notion, arising for the confining potential (5), is the confining dimension 

d d 1 

Z a=l ' 4 a 

The density of states (3), for potential (5), becomes 

Pie) = -^r , (8) 

where T(s) is a gamma-function. Form (8) is valid for any d > 1 and s > 1/2. The parameter 
Id is 

d w l/2+l/n a 

~*r n r(1 a + 1/na) ( 9 ) 

for d > 1. For instance, for d = 1, it is 

7T /CJ XS 



7i = — — ^-^—r~r — (d = 1) ■ (10) 

11 r(i + i/n) V 2 y v ; v ; 

In the case of harmonic traps, when n a = 2, the confining dimension (7) coincides with 
the spatial dimension d = s. Then 7^ = Uq for all d > 1, and the density of states is 
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III. MODIFIED SEMICLASSICAL APPROXIMATION 



Observable quantities, calculated with the density of states (8), are expressed through 
the Bose-Einstein function g s (z), in which s is the confining dimension (7) and z = is 
fugacity, with f3 = 1/T being inverse temperature and /i, chemical potential. Below the 
Bose condensation temperature T c , one has fi — > and z 1. The Bose-Einstein function 
<? s (z) diverges as z — > 1, if s < 1. Those observable quantities that contain g s (l), with 
s < 1, cannot be defined, because of the divergence of g s (l). This imposes the limits of 
applicability for the semiclassical approximation. 

However the divergence of <7 S (1), for s < 1, is related to the usage of the thermodynamic 
limit assuming an infinite system with an infinite number of particles. For this case, integra- 
tions over momenta start from zero, which implies that the minimal momentum is zero. But 
if the system is finite, containing a finite number of particles N, though may be very large, 
and in addition is confined in space by a trapping potential with a characteristic length 
Iq, then the minimal momentum of a particle is not zero, but rather is a finite quantity 
kmin = 1 / Zo- Respectively, the minimal energy is k^ nin /2m, and the related dimensionless 
minimal energy is 

v _PK m _j__uJ 1 nn 

mm 2m ~ 2ml\ ~ 2T " 1 } 

The Bose-Einstein functions g n (z), arising in the process of calculating physical quantities, 
are defined through integrals over the dimensionless energy variable u = f3e. For a finite 
confined system, the integration should start from the minimal value (11), which yields 

g n {z) = — / du. 12 

Y[n) J Umin e u - z 

The lower integration limit, according to inequality (2) is small, 

u min = ^ < 1 , (13) 

though it is not strictly zero. For n > 1, when the integral (12) converges, the value (13) is 
negligible and can be replaced by zero. However, for n < 1, when the integral can diverge 
at z — > 1, the lower integration limit must be kept finite, being given by Eq. (11). In this 
way, we can define the integral (12) for all n, which gives 

2T 

9l (l)^\n— (n=l), 

1 /2T\ 1 ~ n 

9n{l) = -. — (0 < n< 1) , 

(1 - n)T(n) \u J 

2T 

g Q {\) - — (n = 0) . (14) 
For the last equality in Eqs. (14), the relation 

Y z 9n(z) = \ g n -i(z) 
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is used, being valid for all n. For negative values n < 0, the integral (12) has the form of 
the second of Eqs. (14). Note that for all n > 0, the values <? n (l) are positive, while for 
1/2 < n < 0, these values become negative, since then T(n) < 0. 
The total number of particles can be represented as 

N = N + — g s (z) , (15) 
Id 

where N is the number of condensed atoms. From here, the critical temperature T c follows 
as the temperature at which fj, — > 0, z — > 1, and Nq — > 0, which results in 

l/s 

(16) 



T, 



c 



7dN 



M l ). 

For s > 1, one has g s (l) = C( s )j with ((s) being a Riemann zeta function. In the standard 
picture, g s (l) would diverge for all s < 1, which would lead to the conclusion that then 
T c — > 0. That is, then a finite condensation temperature would not exist for one-dimensional 
systems trapped by a potential with the power n > 2, which includes the harmonic potential 
[2,7,11]. But, as is explained above, for a finite confined system, we have to take the value 
g n (l) given by Eqs. (14). 

For s = 1, which happens for a one-dimensional harmonic trap (d — 1, n — 2), we find 
the critical temperature 

which is an immediate consequence of Eqs. (14) and (16). Equation (17) exactly coincides 
with the equation for T c obtained in a purely quantum-mechanical calculation [12]. Iterating 
Eq. (17), with taking account of the inequality 

^ « e 2N (N » 1) , 

UJ 

we get the condensation temperature 

T ^mr } (» = "=i. » = fl («) 

for the one-dimensional harmonic trap. 

In the case of one-dimensional traps, with confining potentials of powers n > 2, we can 
have the confining dimension in the interval 1/2 < s < 1. In such a case, the condensation 
temperature is 

Thus, we come to a conclusion that finite critical temperatures do exist for one-dimensional 
traps for all powers of the confining potential. 

Dealing with harmonic traps, one has s — d, so that for two-and three-dimensional traps, 
the critical temperatures are 

l/d 

(s = d = 2,3). (20) 



T c = u 



N 



((d). 

And for a one-dimensional harmonic trap, T c is given by Eq. (18). 
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IV. EFFECTIVE THERMODYNAMIC LIMIT 



Effective thermodynamic limit for three-dimensional harmonic traps has been defined 
earlier [3]. Here this notion is generalized for all spatial dimensions d > 1 and arbitrary 
confining dimensions s > 1/2. 

The basic idea in defining an effective thermodynamic limit is the common agreement 
that, when the number of particles tends to infinity, N — > oo, the extensive observable quan- 
tities must vary proportionally to N. One of the main extensive observables is the internal 
energy E. Hence, the thermodynamic limit in the most general sense can be understood as 
the limit 

E 

N — > oo , E — > oo , ]y — * const , (21) 
where a nonzero constant is assumed. For the internal energy, we have 

# _ sg s+ i(z) ^ s+1 

N N ld 1 1 

Since s > 1/2, the value of g s+ i(z) is always finite. This means that the limit (21) can be 
rewritten as 

N -> oo , 7 d -> , Nj d -> const . (23) 

This is the most general form of the thermodynamic limit, valid for all power-law potentials 
in any spatial dimensionality. 

Let us consider in more detail the case of unipower potentials, for which n a = n. Then 
the confining dimension (7) becomes 

s=( 1 - + -)d. (24) 



2 n 



The coefficient (9) takes the form 



ld= T^ll/n) {yJ (25) 

where o;o is defined in Eq. (6). This shows that 7 rf oc Uq for all d > 1. Thence, the limit 
(23) can be represented as the effective thermodynamic limit 

N -> oo , cu -> , Nu s -> const . (26) 

The latter implies that 

u oc iV- 1/s (JV -> oo) . (27) 

In a particular case of a three-dimensional harmonic trap, when s = d = 3, the limit (26) 
reduces to the known result [3]. But, generally, the limit (26) is valid for all confining 
dimensions (24). 
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The condensation temperatures, found above, have sense for confined systems with a 
finite, though very large, number of particles N. The notion of the thermodynamic limit 
(26) allows one to analyze the behaviour of T c as N — > oo. To this end, from Eq. (16) we 
have 

T c oc const (s > 1) . (28) 
For the one-dimensional harmonic trap, Eq. (18) gives 

Tc K LTiv (s = 1} • (29) 

And for very low confining dimensions s, from Eq. (19) we find 

T c oc N 1 - 1 / 3 ^ (s < 1) . (30) 

Equations (28) to (30) demonstrate the behaviour of T c as N — > oo. But for finite JV, 
the corresponding values in Eqs. (16), (18), and (19) can be finite and well defined. 



V. SPECIFIC HEAT DISCONTINUITY 

Specific heat, under a fixed number of particles, is given by the derivative 

Cn = not- (31) 

With the internal energy (22), for the temperatures above T c , one has 

C N = s(s+l) 9 -^± - s 2^(£) ( T >T C ). (32) 
9s{z) 9s-i(z) 

One may notice that for a three-dimensional uniform system, when n — > oo and s — > 3/2, 
Eq. (32) reduces to the known expression [13,14]. Below the condensation temperature T c , 
we get 

c ^ a{a+1)i wr(B' {T<T ' ] - (33) 

Expression (33) is finite and positive for all s > 1/2. Defining the specific heat discontinuity 
at the critical point as 

AC N = C^(T C + 0) - C N (T C - 0) , (34) 

we find 

AC N = -s 2 . (35) 

In the standard approach, because of the divergence of (? s _i(l) for s < 2, the specific heat 
discontinuity (35) could be considered [11] only for s > 2. Here we extend the analysis for 
all s. 
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First, for higher confining dimensions, we have 

ACV = s 2 -f^-, (s>2), (36) 

which agrees with the previous result [11]. For lower confining dimensions, which could not 
be treated earlier, we obtain the following expressions: 

2ir 2 

AC » = -3Mzfj^) (s = 2) ' (37) 



2-s 



AC N = -s 2 ((s)(2-s)T(s-l)(^) (Ks<2), (38) 



2T C 

AC « = - m ln O < s = » • < 39 » 

For all s > 1, the specific heat jump is negative, which tells that 

C N (T c -0)>C N (T c + 0) (s>l). 
But the sign of the discontinuity changes for s < 1, demonstrating that 

C N (T c -0) <C N (T c + 0) (s< 1) . 

The jump (40) becomes positive, since F(s — 1) is negative for 1/2 < s < 1. 

In the case of harmonic traps, the specific heat discontinuity ACn is given by Eq. (39) 
for d — 1, by Eq. (37) for d = 2, and for a three-dimensional harmonic trap, one has 

AC N = - - 2 C(3) (s = d = 3) . 
For finite systems, the specific heat is always discontinuous at T c . 

VI. DENSITY FLUCTUATIONS AND COMPRESSIBILITY 

Density fluctuations in any statistical system can be quantified [13,15] by the isothermal 
compressibility 

A 2 (iV) 1 (dp\ ,„, 

where p is a mean particle density, and which is connected with the number-of-particle 
dispersion 



8 



A 2 {N) = <N 2 > - <N> 2 



(42) 



For the latter, one has 



A 2 (iV)=T^ = ^(^) . (43) 




The total number of particles N — N + Aq is a sum of the number of condensed atoms 
and the number of atoms outside the condensate. In a gauge- symmetric grand canonical 
ensemble, the fluctuations of condensate are known to be anomalous [16,17]. However, as 
is also perfectly known [16], this is just an artifact that can be easily removed by breaking 
the gauge symmetry. The most efficient way of gauge-symmetry breaking is by means of 
the so-called Bogolubov shift [18,19], when the field operators of condensed particles are 
replaced by the nonoperator quantities, representing the condensate wave functions. In that 
way, the operator of the number of condensed particles is replaced by its average, as a result 
of which only the noncondensed particles contribute to the dispersion 

A 2 (iV) = A 2 (A\) = T ^ . (44) 

O/jL 

The necessity of introducing broken gauge symmetry in grand canonical ensemble, in order 
to eliminate fictitious condensate fluctuations, was discussed in great detail by ter Haar 
[16] and Hohenberg and Martin [20], and carefully explained by Bogolubov [18,19]. The 
asymptotic exactness of the Bogolubov shift was proven by Ginibre [21]. 
The number of noncondensed particles, according to Eq. (15), is 

N 1 = — g s (z) . (45) 
Id 

From here, the dispersion (44) writes as 

A 2 (Aq) = - g_ x {z) . (46) 
Id 

This, invoking relation (16), can be rewritten as 

A 2 (iV 1 ) = iV^^f-V . (47) 



9s{l) \T C 

Involving the notion of the thermodynamic limit, particle fluctuations can be classified into 
normal and anomalous. This is described more fully in the review article [22]. When 
A 2 (iVi) oc N, fluctuations are called normal, while if A 2 (iVi) oc N a , with a > 1, they are 
termed anomalous. 

At temperatures T > T c , above the condensation point, one has N± — N and z < 1. For 
the confining dimension s > 1, particle fluctuations are normal, since 

A 2 (iVi) oc iV (s > 1, T > T c ) . (48) 

Respectively, the compressibility (41) is finite for all iV — > oo. For lower confining dimensions 
s < 1, we find 
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^ = ( i-,)(i-T)r ( ,-i h , GiT ■ < 49 > 

Then, because T(s — 1) < for 1/2 < s < 1, the compressibility (41) becomes divergent and 
negative, 

k t oc -jV~ 1+1/s (s < 1, T > T c ) . (50) 

Such a behaviour of the compressibility means that the system is unstable. 

Above the condensation temperature T c , the trapped gas is stable only for the confining 
dimensions s > 1. For harmonic traps, for which s — d, the gas is stable in all spatial 
dimensions d > 1. 

The situation is more interesting for the temperatures T < T c below the condensation 
point. Then the dispersion (46) becomes 

A 2 (iV 1 ) = -^i(l). (51) 

Id 

There exists a rich variety of different cases: 

C(*-l)/T 



AW = W -_^_, ( S>2 ), (52) 



TV /T\ 2 , /2T 



A w = «i) Is J ln U> (s - 2) - (53) 



AT /2TA 2_S /T^ 2 



^^»RWi feJ (1<s<2) ' (54) 

A 2 (iV 1 ) = 2(-) 2 ( a = l). (55) 
Wo/ 

And for s < 1, the dispersion A 2 (Ai) has the same form as in Eq. (54). In the thermody- 
namic limit, we find 

A 2 (A\) oc AT (s > 2) , 

A 2 (A r 1 ) oc A^ln A" (s = 2) , 

A 2 (A\) oc N 2 ' s (1< s < 2) , 

A 2 (A\) oc AT 2 ( s = 1) , 

A 2 (A%) oc —N 2 / s (s < 1) . (56) 

Respectively, the behaviour of the compressibility is 

kt oc const (s > 2) , 
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k t oc In N (s = 2) , 
/t T oc N 2 ' 3 - 1 (1< s < 2) , 
kt oc N (s = 1) , 

k t oc -jV 2/8_1 (s < 1) . (57) 

These equations demonstrate that the fluctuations are anomalous for all s < 2. A neg- 
ative compressibility for s < 1 implies strong instability of the system. For the confining 
dimensions in the interval 1 < s < 2, the compressibility is positive but displays nonther- 
modynamic behaviour diverging in the thermodynamic limit. An actual divergence of kt 
happens only for infinite systems, when A" — > oo, which would imply instability. When one 
deals with finite system, with a large, though finite, number of particles N 1, then the 
compressibility does not really diverge, but becomes very large. This means that there exist 
very strong density fluctuations in the system. Such strong fluctuations arise in the case 
of the low confining dimension s < 2, which corresponds to large confining powers n a and 
low spatial dimensions d. Fluctuations are known to be strong in low-dimensional systems 
[23,24]. For higher confining dimensions s > 2, fluctuations are always normal, so that for 

I + S^>2 (58) 

z a=l "« 

the system is stable. 

Since the standard traps are usually harmonic, let us pay a special attention to the 
harmonic confining potentials, when s = d and 7^ = Uq. Then the dispersion (51) becomes 

A 2 (A> 1 )=^ 1 (1)(-)" • (59) 
For different spatial dimensions, we find 

^^"miff ( ^ 3) - (60) 

The last expression agrees with the corresponding result by Politzer [25] for a three- 
dimensional harmonic trap. All other formulas of this section are new. 
Dispersions (60) in the thermodynamic limit behave as 

A 2 (A>0 oc A^ 2 (d = 1) , 

A 2 (A>i) oc N\nN (d = 2) , 
A 2 (A\) oc AT (d = 3) . 
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So that the dispersions for d — 1 and d = 2 are anomalous. The related compressibilities 
possess the limits 

k t oc N (d = 1) , 

kt oc In iV (d = 2), 

k t oc const (d = 3) . 

Anomalous values of the compressibilities for low-dimensional harmonic traps signify insta- 
bility caused by the existence of very strong fluctuations in such traps. 

VII. CONCLUSION 

The semiclassical approximation is generalized, which makes it possible to essentially 
extend the region of its applicability. Bose-Einstein condensation in traps of low confin- 
ing dimensions is described, for which the standard approach could not be used. Trapping 
potentials of arbitrary power laws are considered. Specific-heat discontinuities and isother- 
mal compressibilities are analyzed. Effective thermodynamic limit is defined for any spatial 
dimension and for arbitrary powers of confining potentials. It is shown that the modified 
semiclassical method for harmonic traps yields the results coinciding with those obtained by 
means of quantum-mechanical calculations, when these are available. 
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